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ABSTRACT 

We present results from three-dimensional particle simulations of collisionless shock formation, with rela- 
tivistic counter- streaming ion-electron plasmas. Particles are followed over many skin depths downstream 
of the shock. Open boundaries allow the experiments to be continued for several particle crossing times. The 
experiments confirm the generation of strong magnetic and electric fields by a Weibel-like kinetic streaming in- 
stability, and demonstrate that the electromagnetic fields propagate far downstream of the shock. The magnetic 
fields are predominantly transversal, and are associated with merging ion current channels. The total magnetic 
energy grows as the ion channels merge, and as the magnetic field patterns propagate down stream. The elec- 
tron populations are quickly thermalized, while the ion populations retain distinct bulk speeds in shielded ion 
channels and thermalize much more slowly. The results help reveal processes of importance in collisionless 
shocks, and may help to explain the origin of the magnetic fields responsible for afterglow synchrotron/jitter 
radiation from Gamma-Ray Bursts. 

Subject headings: acceleration of particles, gamma rays: bursts, instabilities, magnetic fields, plasmas 
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1. INTRODUCTION 

The existence of a strong magnetic field in the shocked ex- 
ternal medium is required in order to explain the observed 
radiation in Gamma-Ray Burst afterglows as synchrotron ra- 
diation (e.g. Panaitescu & Kumar 2002). Nearly collisionless 
shocks, with synchrotron-type radiation present, are also com- 
mon in many other astrophysical contexts, such as in super- 
nova shocks, and in jets from active galactic nuclei. At least 
in the context of Gamma-Ray Burst afterglows the observed 
synchrotron radiation requires the presence of a stronger mag- 
netic field than can easily be explained by just compression of 
a magnetic field already present in the external medium. 

Medvedev & Loeb (1999) showed through a linear kinetic 
treatment how a two- stream magnetic instability - a general- 
ization of the Weibel instability (Weibel 1959; Yoon & David- 
son 1987) - can generate a strong magnetic field (e#, defined 
as the ratio of magnetic energy to total kinetic energy, is 10~ 5 - 
10 _1 of equipartition value) in collisionless shock fronts (see 
also discussion in Rossi & Rees 2003). We note in passing 
that this instability is well-known in other plasma physics dis- 
ciplines, e.g. laser-plasma interactions (Califano el al. 1998; 
Yang et al. 1992), and has been applied in the context of pulsar 
winds by Kazimura et al. (1998). 

Using three-dimensional particle-in-cell simulations to 
study relativistic collisionless shocks (where an external 
plasma impacts the shock region with a bulk Lorentz factor 
T = 5 - 10), Frederiksen et al. (2003), Nishikawa et al. (2003), 
and Silva et al. (2003) investigated the generation of mag- 
netic fields by the two- stream instability. In these first stud- 
ies the growth of the transverse scales of the magnetic field 
was limited by the dimensions of the computational domains. 
The durations of the Nishikawa et al. (2003) experiments were 
less than particle travel times through the experiments, while 
Silva et al. (2003) used periodic boundary conditions in the 
direction of streaming. Further, Frederiksen et al. (2003) 
and Nishikawa et al. (2003) used electron-ion (e~p) plasmas, 
while experiments reported upon by Silva et al. (2003) were 
done with e~e + pair plasmas. 
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Here, we report on 3D particle-in-cell simulations of rela- 
tivistically counter- streaming e~p plasmas. Open boundaries 
are used in the streaming direction, and experiment durations 
are several particle crossing times. Our results can help reveal 
the most important processes in collisionless shocks, and help 
explain the observed afterglow synchrotron radiation from 
Gamma-Ray Bursts. We focus on the earliest development 
in shock formation and field generation. Late stages in shock 
formation will be addressed in successive work. 

2. SIMULATIONS 

Experiments were performed using a self-consistent 3D3 V 
electromagnetic particle-in-cell code originally developed for 
simulating reconnection topologies (Hesse et al. 1999), rede- 
veloped by the present authors to obey special relativity and 
to be second order accurate in both space and time. 

The code solves Maxwell's equations for the electromag- 
netic field with continuous sources, with fields and field 
source terms defined on a staggered 3D Yee-lattice (Yee 
1966). The sources in Maxwell's equations are formed by 
weighted averaging of particle data to the field grid, using 
quadratic spline interpolation. Particle velocities and posi- 
tions are defined in continuous (r, 7v)-space, and particles 
obey the relativistic equations of motion. 

The grid size used in the main experiment was (x,y,z) = 
200 x 200 x 800, with 25 particles per cell, for a total of 
8 x 10 8 particles, with ion to electron mass ratio mi/m e = 16. 
To adequately resolve a significant number of electron and 
ion skin-depths (5 e and Si), the box size was chosen such that 
L XJ = 10(5; ~ 40(5 e and L z ~ 40(5/ ~ 1605 e . Varying aspect and 
mass ratios were used in complementary experiments. 

Two counter- streaming - initially quasi-neutral and cold - 
plasma populations are simulated. At the two-stream interface 
(smoothed around z = 80) a plasma (z < 80) streaming in the 
positive z-direction, with a bulk Lorentz factor r = 3, hits an- 
other plasma (z > 80) at rest in our reference frame. The latter 
plasma is denser than the former by a factor of 3. Experiments 
have been run with both initially sharp and initially smooth 
transitions, with essentially the same results. The long sim- 
ulation time gradually allows the shock to converge towards 
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Fig. 1 . — The left hand side panel shows the longitudinal electron current density through a transverse cut at z = 100, with a small inset showing the ion current 
in the same plane. The right hand side panel shows the ion current at z = 600 = 305/, with the small inset now instead showing the electron current. The arrows 
represent the transverse magnetic field. Both panels are from time t = 1200. 



self-consistent jump conditions. Periodic boundaries are im- 
posed in the x- and v-directions, while the boundaries at z = 
and z = 800 are open, with layers absorbing transverse elec- 
tromagnetic waves. Inflow conditions at z = are fixed, with 
incoming particles supplied at a constant rate and with uni- 
form speed. At z = 800 there is free outflow of particles. The 
maximum experiment duration is 480 uj~ l e (where u pe is the 
electron plasma frequency), sufficient for propagating T « 3 
particles 2.8 times through the box. 

3. RESULTS AND DISCUSSIONS 

The extended size and duration of these experiments make 
it possible to follow the two-stream instability through several 
stages of development; first exponential growth, then non- 
linear saturation, followed by pattern growth and downstream 
advection. We identify the mechanisms responsible for these 
stages below. 

3.1. Magnetic Field Generation, Pattern Growth 
and Field Transport 

Encountering the shock front the incoming electrons are 
rapidly (being lighter than the ions) deflected by field fluc- 
tuations growing due to the two- stream instability (Medvedev 
& Loeb 1999). The initial perturbations grow non-linear as 
the deflected electrons collect into first caustic surfaces and 
then current channels (Fig. 1). Both streaming and rest frame 
electrons are deflected, by arguments of symmetry. 

In accordance with Ampere's law the current channels are 
surrounded by approximately cylindrical magnetic fields (il- 
lustrated by arrows in Fig. 1), causing mutual attraction be- 
tween the current channels. The current channels thus merge 
in a race where larger electron channels consume smaller, 
neighboring channels. In this manner, the transverse magnetic 
field grows in strength and scale downstream. This continues 
until the fields grow strong enough to deflect the much heav- 
ier ions into the magnetic voids between the electron chan- 
nels. The ion channels are then subjected to the same growth 
mechanism as the electrons. When ion channels grow suffi- 



Fig. 2. — Electron (top) and ion (bottom) currents, averaged over the x- 
direction, at time t = 1200. 



ciently powerful, they begin to experience Debye shielding by 
the electrons, which by then have been significantly heated by 
scattering on the increasing electromagnetic field structures. 
The two electron populations, initially separated in 7V-space, 
merge to a single population in approximately 20S e (z = 80- 
200) as seen in Fig. 6. The same trend is seen for the ions - 
albeit at a rate slower in proportion to mi/m e . 

The Debye shielding quenches the electron channels, while 
at the same time supporting the ion-channels; the large ran- 
dom velocities of the electron population allow the concen- 
trated ion channels to keep sustaining strong magnetic fields. 
Fig. 1, shows the highly concentrated ion currents, the more 
diffuse - and shielding - electron currents, and the result- 
ing magnetic field. The electron and ion channels are fur- 
ther illustrated in Fig. 2. Note the limited z-extent of the 
electron current channels, while the ion current channels ex- 
tend throughout the length of the box, merging to form larger 
scales downstream. Because of the longitudinal current chan- 
nels the magnetic field is predominantly transversal; we find 
\B z \/\B tot \~lO- l -lO- 2 . 

Figure 3 shows the temporal development of the transverse 
magnetic field scales around z = 250. The power spectra fol- 
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Fig. 3. — Power spectrum of for z = 250 at different times. 
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FIG. 5. — Total magnetic (full drawn) and electric (dashed) energy in the 
box as a function of time. The inset shows a log-log plot of the same data. 




Fig. 4. — Relative electromagnetic energy density e#. The contour color 
plot shows the power in the transverse magnetic field through the box dis- 
tributed on spatial Fourier modes at t = 2400, with the dotted line marking the 
wavenumber with maximum power. Superposed is the spatial distribution of 
eg, averaged across the beam, at t = 2320 (dashed-dotted) and t = 2400 (full 
drawn), highlighting how EM-fields are advected down through the box. 



back after about 10-30 uj~\. It is indeed obvious from the pre- 
ceding discussion that the ion-electron asymmetry is essential 
for the survival of the current channels. 

From the requirement that the total plasma momentum 
should be conserved, the (electro)magnetic field produced by 
the two- stream instability acquires part of the z-momentum 
lost by the two-stream population in the shock; this opens the 
possibility that magnetic field structures created in the shock 
migrate downstream of the shock and thus carry away some 
of the momentum impinging on the shock. 

Our experiments show that this does indeed happen; the 
continuous injection of momentum transports the generated 
field structures downstream at an accelerated advection speed. 
The dragging of field structures through the dense plasma acts 
as to transfer momentum between the in- streaming and the 
shocked plasmas. 



low power-laws, with the largest scales growing with time. 
The dominant scales at these z are of the order Si at early 
times. Later they become comparable to L XJ . Figure 4 cap- 
tures this scaling behavior as a function of depth for t = 2400. 

The time evolutions of the electric and magnetic field ener- 
gies are shown in Fig. 5. Seeded by fluctuations in the fields, 
mass and charge density, the two-stream instability initially 
grows super-linearly it = 80- 100), reflecting approximate ex- 
ponential growth in a small sub-volume. Subsequently the 
total magnetic energy grows more linearly, reflecting essen- 
tially the increasing volume filling factor as the non-linearly 
saturated magnetic field structures are advected downstream. 

At ttt 1 100 the slope drops off, due to advection of the gen- 
erated fields out of the box. The continued slow growth, for 
t > 1 100, reflects the increase of the pattern size with time (cf. 
Fig. 3). A larger pattern size corresponds to, on the average, 
a larger mean magnetic energy, since the total electric current 
is split up into fewer but stronger ion current channels. The 
magnetic energy scales with the square of the electric current, 
which in turn grows in inverse proportion to the number of 
current channels. The net effect is that the mean magnetic 
energy increases accordingly. 

The magnetic energy density keeps growing throughout our 
experiment, even though the duration of the experiment (480 
significantly exceeds the particle crossing time, and also 
exceeds the advection time of the magnetic field structures 
through the box. This is in contrast to the results reported by 
Silva et al. (2003), where the magnetic energy density drops 



3.2. Thermalization and Plasma Heating 

At late times the entering electrons are effectively scattered 
and thermalized: The magnetic field isotropizes the velocity 
distribution whereas the electric field generated by the e~-p 
charge separation acts to thermalize the populations. Figure 
6 shows that this happens over the ~ 20 electron skin depths 
from around z = 80 - 200. The ions are expected to also ther- 
malize, given sufficient space and time. This fact leaves the 
massive ion bulk momentum constituting a vast energy reser- 
voir for further electron heating and acceleration. Also seen 
in Fig. 6, the ions beams stay clearly separated in phase space, 
and are only slowly broadened (and heated). 

We do not see indications of a super-thermal tail in the 
heated electron distributions, and there is thus no sign of sec- 
ond order Fermi- acceleration in the experiment presented in 
this Letter. Nishikawa et al. (2003) and Silva et al. (2003) re- 
ported acceleration of particles in experiments similar to the 
current experiment, except for more limited sizes and dura- 
tions, and the use of an e~e + plasma (Silva et al. 2003). On 
closer examination of the published results it appears that 
there is no actual disagreement regarding the absence of ac- 
celerated particles. Whence, Nishikawa et al. (2003) refer to 
transversal velocities of the order of 0.2c (their Fig. 3b), at a 
time where our experiment shows similar transversal veloci- 
ties (cf. Fig. 6) that later develop a purely thermal spectrum. 
Silva et al. (2003) refer to transversal velocity amplitudes up 
to about 0.8c (their Fig. 4), or ~ 2, with a shape of the 
distribution function that appears to be compatible with ther- 
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Fig. 6. — Thermalization and longitudinal acceleration, illustrated by scat- 
ter plots of the electron (orange) and ion (blue) populations. Note the back- 
scattered electron population (v z j(v) < 0). 



mal. In comparison, the electron distribution illustrated by 
the scatter plot in Fig. 6 covers a similar interval of vj, with 
distribution functions that are close to (Lorentz-boosted) rel- 
ativistic Maxwellians. Thus, there is so far no compelling 
evidence for non-thermal particle acceleration in experiments 
with no imposed external magnetic field. Thermalization is a 
more likely cause of the increases in transversal velocities. 

Frederiksen et al. (2003) reported evidence for particle ac- 
celeration, with electron gammas up to ~ 100, in experi- 
ments with an external magnetic field present in the up-stream 
plasma. This is indeed a more promising scenario for parti- 
cle acceleration experiments (although in the experiments by 
Nishikawa et al. 2003, results with an external magnetic field 
were similar to those without). Figure 6 shows the presence of 
a population of back- scattered electrons (v z j < 0). In the pres- 
ence of an external magnetic field in the in- streaming plasma, 
this possibly facilitates Fermi acceleration in the shock. 

4. CONCLUSIONS 

The experiment reported upon here illustrates a number of 
fundamental properties of relativistic, collisionless shocks: 
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1. Even in the absence of a magnetic field in the up- 
stream plasma, a small scale, fluctuating, and predominantly 
transversal magnetic field is unavoidably generated by a two- 
stream instability reminiscent of the Weibel-instability. In the 
current experiment the magnetic energy density reaches a few 
percent of the energy density of the in-coming beam. 

2. In the case of an e~p plasma the electrons are rapidly 
thermalized, while the ions form current channels that are the 
sources of deeply penetrating magnetic field structures. The 
channels merge in the downstream direction, with a corre- 
sponding increase of the average magnetic energy with shock 
depth. This is expected to continue as long as a surplus of bulk 
relative momentum remains in the counter- streaming plas- 
mas. 

3. The generated magnetic field patterns are advected 
downstream at speeds intermediate of the streaming and rest- 
frame plasmas. The electromagnetic field structures thus pro- 
vide scattering centers that interact with both the fast, in- 
coming plasma, and with the plasma that is initially at rest. As 
a result the electron populations of both components quickly 
thermalize and form a single, Lorentz-boosted thermal elec- 
tron population. The two ion populations merge much more 
slowly, with only gradually increasing ion temperatures. 

4. The observed strong turbulence in the field structures 
at the shocked streaming interface provides a promising envi- 
ronment for particle acceleration. 

We emphasize that quantification of the interdependence 
and development of ejj and eg is accessible by means of such 
experiments as reported upon here. 

Rather than devising abstract scalar parameters e# and €u, 
that may be expected to depend on shock depth, media den- 
sities etc., a better approach is to compute synthetic radiation 
spectra directly from the models, and then apply scaling laws 
to predict what would be observed from corresponding, real 
supernova remnants and Gamma-Ray Burst afterglow shocks. 
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